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Abstract 

We present a consistent analysis of T sum rules and i?-meson semileptonic width in 
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1 Introduction. 



The rich phenomenology of the heavy quark mesons has been recognized as a source of 
information on strong interaction dynamics and a clean place to determine fundamental 
parameters of the standard model such as the strong coupling constant, heavy quark masses 
and mixing angles. The necessity of the quantitative analysis of the strong interaction effects 
in the heavy flavor physics leads to developing new theoretical methods such as nonrelativistic 
QCD (NRQCD) [||] for description of heavy quarkonium and heavy quark effective theory 
(HQET) for description of heavy-light mesons. One of the most promising applications 
of these techniques to the bottom quark physics is the analysis of the T sum rules and 
5-meson semileptonic decays (see |^ for a recent review). The strict estimate of the 
actual precision of theoretical predictions is now becoming important since the experimental 
precision is rather high and is quantitatively comparable with the last available terms of 
perturbative approximation. Recently a significant progress has been achieved in calculations 
of the high order corrections for the T sum rules |, ||, [1^ and inclusive semileptonic 

width of S-mesons [|ll|, [T|, |l| . 

In this paper we present a simultaneous analysis of T sum rules and the S-meson in- 
clusive semileptonic width Fsi in the next-to- next-to-leading order (NNLO). On the basis 
of this analysis we determine the mass parameter of the bottom quark in the second order 
of perturbation theory^ and give a new accurate estimate of Cabibbo-Kobayashi-Maskawa 

(CKM) matrix element \Vcb\- 

^ Because quark states have not been observed as asymptotic free particle states the notion of their mass 
can not be unambiguously defined. There are different ways to reflect the fact the 6-quark is heavy by 
introducing its mass fl^ . In the present paper we use the notion of the 6-quark pole mass defined strictly 
within the finite order perturbation theory. Within this definition this mass parameter is equivalent (and 
can be reexpressed through a finite series in a^) to any other perturbative definition of the quark mass 
parameter, in particularly, MS mass in every order of the perturbation theory as a parameter of perturbative 
QCD. 
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The results of the analysis of the T system in the NNLO approximation in the coupling 
constant and nonrelativistic expansion have been briefly reported in ref. 0. Then two more 
papers with NNLO analysis of the T system appeared |l^ . Since there is some difference 
between the analyses of ref. § and refs. in some details of the treatment of the 

problem, in this paper we give more extensive discussion paying the main attention to the 
self-consistence of the approach, uncertainties of theoretical estimates and the reliability of 
results. We discuss also the general structure and the asymptotic character of perturbative 
expansion for the moments of the T system spectral density. By studying the behavior of the 
expansion around the Coulomb approximation for different moments we estimate the critical 
order where the asymptotic growth of the coefficients of the series in as for the moments 
starts. 

The obtained moments of the T system spectral density are related then to the inclu- 
sive 5-meson semileptonic width up to the NNLO order. In this way we avoid the strong 
dependence of Fgi on m;, and reduce the theoretical uncertainty in |Vcf,|. In particular, the 
obtained perturbative expansion for Vcb converges well (in heuristic sense) up to NNLO in 
as. 

The paper is organized as follows. In the next Section we describe the vacuum polarization 
function of heavy quarks near threshold in the NNLO. In Section |^ we discuss the T sum rules 
and numerical estimates for the 6-quark pole and MS mass within perturbation theory. In 
Section^ the inclusive i?-meson semileptonic width and \ Vcb\ matrix element are determined. 
The last Section contains our conclusions. In Appendix we give an explicit analytical formula 
for the polarization function of heavy quarks near threshold. 
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2 Next-to-next-to-leading order heavy quark vacuum 
polarization function near threshold. 

We study the near threshold behavior of the vacuum polarization function n(s) of the 6-quark 
vector current = 67^6 

{q,q. - 9,uq') n(g2) = ^ J da:e^''-(0|Tj^(a;)j,(0)|0) . 

The finite order perturbation theory expansion in the strong coupling constant as is known 
to break down near the two-particle production threshold, s ~ Ami where rrib is the 6-quark 



pole mass. However, near the threshold, i.e. for a small 6-quark velocity v = \Jl — Aml/s <ti 
1, the nonrelativistic approximation becomes valid DTf. In the NNLO the nonrelativistic 
expansion for the polarization function 11 (s) has the form 

N / A k'^ \ 

n(s) = TT^ C,(a,)G(0, 0, k) + -—Gc{0, 0, k) (1) 

with k = \Jm1 — s/A being a natural energy variable near threshold. The quantity Ch{oLs) is 
a perturbative coefficient that matches correlators of relativistic and nonrelativistic vector 
currents and accounts for the hard part of the QCD corrections. It is now known to the 
second order in as expansion 

Cn{as) = l-ClCp'^ + ClcJ'^'' 



with = 4 ig and 



Cl=\ C 3 + n 2 \Cf-\ + — C 3 + n 2 Ca 

^ \A ^^^ 3 18/ V 36 2 ^^^ 3 72 J 

where as is defined in MS renormalization scheme with the scale parameter /x. Another 
scale parameter /ij is a factorization scale which separates contributions coming from the 



hard and soft momentum regions and plays the role of an infrared cutoff in the diagrams 
contributing to the quantity Ch{as) 0- The color symmetry SU{3) group invariants 



for QCD are Ca = 3, Cp = 4/3, Tp = 1/2, nj is the number of light fermion flavors, and 
Po = llCyi/3 — 4Tpnj/3 is the first /5-function coefficient. Here •je = 0.577216 ... is the Euler 
constant and ({z) is the Riemann ^-function. The quantity G(x, y,A;) is the nonrelativistic 
Green function. This Green function sums up the singular threshold corrections and satisfies 
the following Schrodinger equation 



Ax A, 



rUb 



4m? 47r KiirJ 



+AnaV{x) + AbfV{^, d^, S) + — ) G(x, y, k) = 5(x - y) (3) 

where Ax = d^, Vc{x) = —Cpds/x is the Coulomb potential, x = |x|, AnaV{x) = 
—CACF0i1/{2mbx'^) is the so called non-Abelian potential of quark-antiquark interaction [|18|, 
AbfV{^, (9x, S) is the standard Breit-Fermi potential (up to the color factor Cp) containing 



the quark spin operator S, e.g. [jT9|. The terms AiV {i = 1, 2) represent the first and second 



order perturbative QCD corrections to the Coulomb potential ||2^, |2T|, [22| pl 



A,V{x) = ^Vc{x){C', + Cl In(x^)), 
47r 



^2V{x) = (^^j Vc{x){C', + C^Hxf,) + Ciln'ixf,)) (4) 



where 



=ai + 2(5olE, Cl = 2A 



0, 



^2 = ( y + 47I j f3^, + 2(A + 2f3oaihE + as, 



^The value of the 02 coefficient in cq. obtained in rcf. [£1| exceeds the correct resuh obtained recently 
p2t by 27r^C^. Though the value of ref. |2l] was used in the previous analysis of the T sum rules and 
B-meson semileptonic width its difference from the present one is numerically small and results in no 
significant change for the numerical estimates of nii,, ag and \ Vcb\ given in these papers. In the present paper 
the most recent value of ref. [|2| is used. 
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31^ 20 

/4343 . TT^ 22 ^2 /1798 56 



V 162 

y - 16C(3)) CpTprif + (yr^n/)' , 

The second term in eq. (|ID is generated by the operator of dimension five in the nonrelativistic 
expansion of the vector current (see, e.g. [^). It contains the Green function of the pure 



Coulomb Schrodinger equation at the origin Gci'^,y, k)\x,y=o In the short distance 

limit X the Coulomb Green function Gc{'^,0, k) has 1/x and Inx divergent terms. 
These terms, however, have no imaginary part and do not contribute to the spectral density 
of the polarization function. Hence they can be subtracted. After the subtraction the 
(renormalized) Coulomb Green function takes the form 

auo, 0. (±) ^)) ,5) 

where \E'„(-2) = d"'\nr{z)/dz^ and T{z) is the Euler F-function. 

The solution of eq. (|^) can be found within the standard nonrelativistic perturbation 
theory around the Coulomb Green function as a leading order approximation 

G(0, 0, k) = Gc(0, 0, k) + AG'(0, 0, k), 
A2 . a. 



AG{0, 0,k) = -J Gc(0, X, k) (^-^ + ^AiV{x) + Gc(x, 0, k)d^ +... (6) 
= Aa^na,bfG + AiG + Af G + A^'^G + . . . 



The corrections to the Coulomb Green function at the origin due to A^, V^a and Vbf 
terms are known analytically |jT6|, The nontrivial part of the calculation consists in 



a proper treatment of some integrals in eq. (^) that correspond to these corrections and 

6 



diverge, or become ill-defined at small x. This divergence is a consequence of the fact that 
the nonrelativistic approximation is not relevant for the description of the short distance 
effects. Moreover, in contrast to the next-to-leading order there exists a divergence in the 
imaginary part of the Green function in NNLO that contributes to the spectral density. The 
divergence can be regularized by introducing an ultraviolet cutoff fif. Then, following the 
general line of the effective field theory approach one has to match the calculation of the 
Green function to the calculation of the coefficient Chicks) i-^- to match the regularization 
procedures for the G(0, 0, k) and Ch- This can be done by comparing eq. (|T]) and the result 
of perturbative calculation of the spectral density [^] in the formal limit <^ f ^ 1 up to 
the order a1 ||I6|, [l^. The result readsQ 

47r ml \ \^if ) 

The next-to-leading (NLO) correction AiG in eq. (|^) due to the first iteration of AiV term 
of the QCD potential and the NNLO correction A2 G due to 1^2^ part of the potential 
have been found in ref. 0]. The correction A2^''G due to the second iteration of AiV term 
which has to be kept in NNLO approximation has been obtained in ref. . We describe the 
details of this calculation below. 

We use the following partial wave representation for the Coulomb Green function 

00 

G'c;(x,y, A;) = Y.{21 + l)Gi{x,y,k)Pi{{^y)/xy) 
1=0 



•^Calculating the NNLO corrections to the nonrelativistic Green function and the hard renormalization 
coefRcient one encounters the divergences with the specific form depending on the regularization procedure. 
They are not included to eqs. (|^, |^ because the divergent terms in the hard renormalization coefRcient cancel 
the divergence in the corrections to the nonrelativistic Green function leaving the logarithmic dependence of 
eq. (H) on the cutoff /U / which compensates corresponding dependence of eq. (Q) . 
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2tt ' ' ' ^o(m + / + l-a,CFm6/(2fc))(m + 2/ + l)! 

where Pi{z) is a Legendre polynomial and L'^{z) is a Laguerre polynomial 

Only / = component of eq. (|^) is necessary for the calculation of the corrections to the 
Green function at the origin up to NNLO approximation. 

Let us consider first the NLO correction. It can be written as 

AiG = ^ (^] £ HMHin) J e-''^Ll{2kx)Ll{2kx)^,V{x)d^ (9) 



m,ri=0 



where 



and we use the equality 



H (ml = m + 1 - 



2k 

Tim + a + 11 



r(a + i)r(m + !)■ 

The integrals in eq. @ corresponding to /S.iV term also diverge at small x. The divergent 
part however is k independent and does not contribute to the spectral density. As a conse- 
quence no matching is necessary for calculation of these corrections. In the representation (^) 
the divergence of the integral at small x is transformed to the divergence of the sum. The 
divergent part of eq. can be separated by the following method. Eq. can be rewritten 
in the form 

2 



^0 J {2kxy J 



where 



r [m) = — — m + 1 — 



m + l)2A;V 2k 
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and we used the property of the Laguerre polynomial 



E 

m=0 



m + a 



In eq. (|10D all sums are convergent and the divergence is contained in the last term. Two 
divergent integrals in the last term of this equation after a regularization take the form 



/ 



e-^'=-ln(/.x)^^ = -lEm + h{kf + ..., 



X 



-2kx 



X 



-dx = L{k) + 



where L{k) = —\n{2k/fi) and ellipsis stands for inessential k independent divergent parts. 
Two finite integrals in eq. (P) are explicitly given by the following expression 



/ e 'Lliz)Ll{z) \n{z)zdz 



(m + l)^i(m + 2), m = n, 
n + 1 



m > n. 



m — n 



I e-'Ll{z) \Yi{z)dz = - ^i(m + 1). (11) 

To compute the first integral we rewrite it in the form containing a derivative with respect 
to an auxiliary parameter e 

J e-^Ll{z)Ll{z)Hz)zdz = ^ (| e'^ Ll{z)Ll{z)z^+^ dz 

By using the relations 



(12) 



e=0 



^ r(/3-a + n) 
„^or(/3-a)r(n + l) — ^ ^' 

V{m + a + l) 



e-'mz)Ll{z)z''dz = 5.. 



r(m+ 1) 



for /3 = l,a = l + £: the integration in the right hand side of eq. (12) can be performed 
analytically. Then taking the derivative in e at e = we get the first line of eq. ([iTl) . The 
second integral in eq. ([TT| ) can be computed using the same technique. Thus the final result 
for the correction is 



AiG = (f^ F\m){m + 1) {C, + {L{k) + ^,{m + 2))Cl 

\m=0 



47r 47r 



oo m-1 



-2 ^ ^ F{m)F{n)——C\ + 2 ^ F(m) (^o^ + {L{k) - - + l))C\) 

m=l n=0 ^ ^ m=0 

+mCl + (-7i.i^(A:) + \L{kf) Ct) . (13) 

Note that the nontrivial part of the calculation is to find the correction to the Green func- 
tion due to the logarithmic correction to the Coulomb potential (^. Computation of the 
correction due to the constant part of the correction to the Coulomb potential is trivial. It 
can be also found from the leading order result by changing the parameter of pure Coulomb 
solution as o.s{l + Clag/'^'n). 

The A'^^G correction has been obtained in the same way. The aS^^G part is finite and 
requires no regularization. It can be computed directly using the representation (H) and the 



first integral of eq. ([TT|) . The results of the calculations are given in Appendix. 

The Green function at the origin can be written in the form which includes only single 
poles in the energy variable. Such a form looks more natural for the Green function of a 
nonrelativistic Schrodinger equation 

where E = —k'^/mb, ipm.E'iO) is the wave function at the origin, the sum goes over bound 
states and the integral is over the state of continuous part of the spectrum. In this way the 
corrections to the Green function stemming from the discrete part of the spectrum reduce 
to corrections to the Coulomb bound state energy levels 



En 



- \^r, (l + AiEm + A^2M A BpEm + Ag^-'E'm + A''2'' Em) (15) 

4(m + 1)^ ^ ^ 



and to the values of Coulomb bound state wave functions at the origin 



l^"^(0)P = 77^^ (1 + Ai^^ + A,5^^ + A^.,^^,^^^^ + A^V^ + A«^^) (16) 



?7r(m + 1)^ 

where Ad5ipm is the correction due to the second term in eq. (|l]). 
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In NLO an explicit analytical expression for the corrections to the bound state parameters 
has the form 

AiE™ = ^ (2^0^ + 2(L(m) + ^^im + 2))C|) , (17) 

(y / 2 

= / + (3L(m) - 1 - 27s + + ^i(m + 2) - 2(m + l)^2(m + l))^!^ 



where 



L m = In — 



The expressions of the NNLO corrections are rather cumbersome and given in Appendix. 

The principal possibility of relating results of eqs. (|1^, |16]) to the mass and leptonic width 
of corresponding T resonances with an account for nonperturbative corrections was discussed 
in H. 

Thus we have described some details of obtaining the complete analytical expressions for 
the vacuum polarization function of heavy quarks near the two-particle threshold in NNLO 
presented in [|^, |^ Q Note that our result does not determine the additive renormalization 
constant in the real part of the polarization function which has to be fixed according to 
the standard normalization condition 11(0) = 0. This constant, however, is inessential for 
physical applications. 



3 T sum rules and 6-quark mass 

The result for the near threshold behavior of the vacuum polarization function of the heavy 
quark is now applied to the analysis of sum rules for the T system. Within the sum rules 
approach the moments A^^'* 

R{s)ds 



(4m|)" 



s=0 



"'The NNLO vacuum polarization function of heavy quarks was also recently obtained in a different 
representation in ref. where the "low scale" mass of 5-quark was studied in the context of T sum rules. 
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of the spectral density R{s) = 127rlmn(s+ie) of the theoretical vacuum polarization function 
of the heavy quarks n(s) should be compared with the experimental ones 

j^exp ^ (4mg)" p Rbis)ds 

Ql Jo ^ ^ 

under the assumption of quark-hadron duality. Here = — 1/3 is the 6-quark electric charge 
and the theoretical parameter rrib is included for convenience. The experimental moments 
j^exp generated by the function Rb{s) which is the normalized cross section of e"*"e~ 
annihilation Rb{s) = a{e~^e~ — hadronS{,5)/o"(e"'"e~ /^^/W") and, at least in principle, can 
be directly found from experiment for any s. In practice, however, the spectral density is 
well measured from experiments only for values of energy rather close to threshold where 
the well pronounced T resonances exist. Therefore numerical values for the experimental 
moments are obtained by saturating the integral in eq. (|18D with the contribution of the 
first six T resonances while the tail of the spectral density at large energy is approximated 
theoretically 

The leptonic widths Tk and masses Mk (A; = 1 ... 6) of the resonances are known with good 
accuracy |2^. For example, for large n the dominant contribution to the moments comes from 



the first T resonance for which rT{i5) = 1.32 ±0.05 keV and Mr{is) = 19460.37 ±0.21 MeV. 
The electromagnetic coupling constant is renormalized to the energy of order nih with the 
result oi'Q^j^{mb) = 1.07a^. The rest of the spectrum beyond the resonance region for energies 
larger than sq ~ (11.2 GeV)^ (continuum contribution) is approximated by the theoretical 
spectral density multiplied by the parameter 0.5 < t < 1.5 which accounts for the uncertainty 
in the experimental data in this energy region. The allowed range for the parameter t is 
chosen to be sufficiently large that gives a very conservative estimate of uncertainty caused 
by all other contributions but the resonance one. For instance, the numerical change of the 
experimental moments produced by a manifest account for the physical continuous spectrum 
due to the open BB production at s > (10.6 GeV)^ is well within the error bars introduced 
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by variation of the t parameter. 

The ordinary perturbative expansion in the couphng constant taken up to a finite order 
is not accurate enough for sufficiently high moments since the coefficients of the expansion 
in low orders of perturbation theory grow fast with the order of the moment. Therefore 
a kind of resummation of the terms giving a dominant contribution is required. These 
terms are identified as those connected with a strong Coulomb interaction in the final state 
of quark-antiquark production. The importance of the Coulomb resummation for large n 
can by directly inferred from the explicit expression for the moments in a pure Coulomb 
approximation 



Ur(i)r(,! + i) n 2 ; r(i=fi)r(n + i; I * ' 



The terms of the series in eq. ( pO]) first increase numerically and then decrease. The number 
m of the term with a maximal numerical magnitude is determined from the relation 

In eq. (^) we neglect a weak m dependence of ({m) for large m. For several first terms with 
m < n and large n the coefficients of the series increase by ^/n in subsequent orders of the 
expansion that makes the Coulomb resummation necessary. Note that for large m > n and 
fixed n the character of asymptotic behavior in m changes and terms become proportional 
to 

(^)"c(")(^)" m 

and the series ( pOD is convergent in at CpOis/'^, < 1. 

On the other hand, no resummation is necessary for small n and finite order perturbation 
theory is sufficient for the phenomenological analysis giving a reasonable precision. Moreover, 
the nonrelativistic approximation is not relevant for the analysis of low moments. The 
ordinary perturbation theory expressions for several first moments of the spectral density 
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are now available with accuracy [0. These low moments, however, cannot be used in 
theoretical formulas for sum rules directly because they get a sizable contribution from the 
large momentum region far from threshold. In this case, the experimental moment necessary 
for comparison can not be found with sufficiently high precision because the spectrum is not 
well known experimentally. 

For numerical estimates it is important to fix the allowed range for the normalization 
point which is present in the explicit formula of the polarization function. Following the 
general line of the renormalization group approach the normalization point has to be chosen 
to minimize the higher order corrections. In fact, the normalization points of as entering 
the coefficient Ch and the nonrelativistic Green function can be different when NNLO cor- 
rections are considered. The difference between the normalization points of the hard and 
soft corrections can be noticed only in higher orders of perturbative expansion. At first sight 
this gives an additional possibility to improve the convergence of the perturbation theory for 
the moments. The typical hard scale of the problem is the heavy quark mass m;,. Indeed, 
one can see that for ^ ~ the NNLO correction to Ch is small in comparison with the 
NLO one. However a naive estimate of the soft "physical scale" as a characteristic scale of 
the Coulomb problem yU ~ mi^as 0, |^ is not acceptable since the direct calculation of the 
NNLO corrections shows that the perturbation theory series for the moments blow up for 
this value of the soft normalization point in MS renormalization scheme. This phenomenon 
can be clearly seen from Fig. 1 and Fig. 2 where the relative weights of the NLO and NNLO 
corrections to the parameters of the nonrelativistic Green function are plotted. Though the 
NLO corrections reach their minimal magnitude at /i ~ mi,as the NNLO corrections to the 
Green function and the moments are completely out of control at this point. At first glance 
this seems to contradict our physical intuition. However since the normalization scale is 
defined in the rather artificial MS scheme the connection of which with hh physics is not 
straightforward there is no reason for coincidence of parameter with any physical scale of 
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the process. The relative weight of the NNLO corrections is stabihzed at /i ~ mf,. Moreover, 
here the fi dependence of the moments is minimal so in our opinion there is no reason to 
split the hard and soft normalization points and use the smaller numerical value for the 
soft normalization scale. The dependence of the moments on yUj is rather weak and we put 

/X/ =/i. 

The range of n which can be used for reliable estimates is also quite restricted. Indeed, 
the low moments cannot be used in sum rules because of the large uncertainty on the 
experimental side due to the poor knowledge of the spectral density at large energies as has 
been already pointed out. From Table |^ one sees that for n^8 the experimental moments are 
rather sensitive to the form of the continuos spectrum beyond the resonance region. Note 
that we assume rather large uncertainty of the continuum to be on the safe side. However, its 
contribution is essentially suppressed in comparison with the resonance one and the resulting 
error of the whole quantity in eq. (|19D is of the same order of magnitude as the uncertainties 
introduced by the resonance contribution. 



t 




M17 


MIT 


0.5 


1.185 


1.015 


1.057 


1.0 


1.315 


1.028 


1.059 


1.5 


1.446 


1.041 


1.061 



Table 1: Sensitivity of the experimental moments to the continuum contribution above Sq 
for So = (11.2 GeV)2. 

On the other hand, the leading nonperturbative power correction to the polarization func- 
tion due to gluonic condensate which can introduce a large uncertainty on the theoretical side 
of sum rules is known to be important for n > 20 The smallness of the nonperturbative 
contribution happened to be the only practical restriction on maximal allowed n in NLO 
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analysis In NNLO analysis, however, the upper hmit on n is stronger and is connected 
with the behavior of the perturbative expansion for the moments as well. From Tables ^, ^ 
one sees that for n>jl2 the perturbative series for the moments is not well convergent and 
the /i dependence of the moments becomes strong that can be considered as an indication 
that the higher order corrections are large here. Perhaps this can mean that the Coulomb 
solution is not a good leading order approximation for these n. A possible explanation of 
this phenomenon can be obtained from the analysis of an asymptotic character of the series 
in Us for the moments. Indeed, the Coulomb resummation (^) extracts from the coefficients 
of perturbation theory for the moments the part that forms a convergent series. This means 
that resummation cannot change the analytic structure of the moment as a function of the 
coupling constant. On the other hand the Coulomb approximation for the rath moment is 
saturated by the large terms of the order m ~ m where m is the solution of eq. (p^. Now 
suppose that the full perturbative expansion for the moments is an asymptotic series with 
mopt being the critical order where the series starts to diverge. In this case if the number 
of the moment is large enough i.e. if m is close to rriopt the asymptotic growth of the high 
order coefficients of the full series becomes more important than the resummed Coulomb 
part of the series. For such n the Coulomb approximation does not saturate the full series 
coefficients in the orders which are dominant in the Coulomb approximation itself and large 
corrections naturally appear. The condition for applicability of the Coulomb solution as a 
leading order approximation now is m < mopt- Taking into account that the perturbative 
expansion around the Coulomb solution is not well convergent for n^l2 from the above 
condition and eq. (^if ) for the physical value of as in MS scheme one finds niopt ~ 3 — 4 
which is a reasonable value for a QCD series. 

Here a remark concerning a specific form of Green function is in order. The difference be- 
tween two representations (|]) and for the Green function appears only in higher orders in 
the coupling constant. We do not attemt to compute the spectral density at positive energies 
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Approx. 


Mf 






Mf, 


LO 


0.9319 


0.6679 


0.5602 


0.5169 


NLO 


0.7905 


0.6633 


0.6413 


0.6682 


NLO* 


0.8123 


0.7235 


0.7590 


0.8694 


NNLO 


1.047 


0.973 


1.026 


1.150 


NNLO* 


1.147 


1.297 


1.760 


2.601 



Table 2: The 0th, 1st and 2nd order theoretical moments for ^^(M^) = 0.118, rrif, = 4.8 GeV 
and fi = ruh. The continuum contribution above sq is subtracted. Here the star (*) stands 
for the moments obtained by using the representation (p^) for Green function. 



in the threshold region because the nonperturbative effects are large and prevent point-wise 
evaluation of this quantity. We are interested in the polarization function in Euclidean re- 
gion where some high order derivatives can be obtained. The nonrelativistic Green function 
has no immediate meaning for us as a source for the spectral density at positive energies 
in threshold region and only serves as a tool for resummation of a special contributions of 
perturbation theory expansion. Therefore these two representations are equivalent in NNLO. 



We find, however, that when the representation (|Tj) is used the perturbative series for the 
moments are more divergent especially for large n (see Table |^) and the dependence of the 
results on n and fi is stronger. Therefore we consider the representation as a preferable 
one. The fact that resummation of the corrections keeping the Coulomb form of the nonrela- 
tivistic Green function spoils the properties of the perturbative expansion can be considered 
as another indication that for large n the Coulomb resummation is not enough and a more 
adequate leading order approximation is necessary with a proper large n behavior. 

Now we turn to a discussion of numerical values of parameters extracted from the analysis 
of sum rules. The sum rules for the T system are not very sensitive to as so NLO [0 and 
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fi (GeV) 


Mf 






KAth 
-'^'20 


6 


0.870 


0.764 


0.769 


0.827 


rrib 


1.047 


0.973 


1.026 


1.150 


4 


1.259 


1.220 


1.334 


1.543 


3 


1.761 


1.867 


2.209 


2.735 



Table 3: The scale dependence of the theoretical moments in the NNLO approximation for 
as{Mz) = 0.118 and nib = 4.8 GeV (the continuum contribution above Sq is subtracted). 



NNLO p analyses give a rather rough estimate 

asiMz) =0.118 ±0.006 



(23) 



which is in agreement with other available data [|28| . This result is obtained by simultaneous 
fit for and m^,. The central values are found by the standard least method [Q. The 
uncertainty of the fit is estimated by fixing to its central value and varying fi and n within 
allowed range. On the other hand, the sum rules are much more sensitive to the 6-quark 



mass so it is instructive to fix ag to the "world average" value as{Mz) = 0.118 |2^ and then 
extract rrib- The final estimate of the bottom quark pole mass is0 



TTib = 4.80 ± 0.06 GeV. 



(24) 



The uncertainty corresponds to the interval 3.5 GeV < /i < 6.5 GeV for 8 < n < 12. In 

fact the central value changes slightly in a wider interval 5 < n < 20 but we restrict the 

range of n to minimize the uncertainty related to fi and t dependence. Note that the scale 

dependence of the moments is mainly due to the scale dependence of as(/i) while the explicit 

dependence on /i is rather weak. 

^In ref. |^ the interval 4.74 GeV < mb < 4.82 GeV has been obtained for 6-quark pole mass. In the present 
paper we use the different range of n and /i for the analysis of the sum rules and give more conservative 
estimate of the uncertainty. 



The perturbation theory relation between the perturbative pole mass and the mass mf,(/i) 
defined in the MS renormalization scheme is known up to two- loop level It is therefore 
straightforward to find the numerical value for the MS mass of 6-quark using result of eq. ( ^if ) 

mb(mfc) =4.21 ±0.11 GeV 

where we assume the interval (|23|) for the strong coupling constant. An interesting fact is 
that the convergence of the perturbative series for MS mass obtained order by order from 
the T sum rules 

mb{mb) = (rnb(mb))^°(l - 0.085 - 0.021) 
is much better than for the pole mass 

mfe = m^^(l- 0.001 + 0.021) 

though the second order terms in the above expansions are rather close numerically. 

Let us emphasize that the convergence of the perturbation theory for the vacuum polar- 
ization function of heavy quark near the threshold and to the moments is not fast in the 
MS renormalization scheme. We have found the NNLO corrections to exceed the NLO ones. 
Furthermore, in the case of 6-quark the corrections due to the perturbative modification of 
the Coulomb instantaneous potential {i.e. related to AGi and AGq^ terms) dominate the 
total correction in the NLO and NNLO. Inclusion of these corrections is quite important for 
consistent analysis of sum rules for the T system. A conjecture that this fact is a conse- 
quence of the asymptotic character of the series which leads to the intrinsic ambiguity in 



the heavy quark pole mass was studied in detail in the literature [|10], ^ ^ It mainly 



based on consideration of the renormalon contribution pO| , |3^ . Within this picture the high 
order contributions to the moments are saturated by the corrections that are associated 
with (a^/^o ln(/ir))'" terms in perturbative series for the heavy quark potential. These cor- 
rections are generated by terms stemming from the running of the coupling constant that 
can be found from renormalization group analysis. However, the explicit calculation shows 
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that in NNLO the corresponding contribution (the one proportional to the coefficient C|) 
provides only ~ 10% of the total NNLO correction i.e. the NNLO correction turns out to 
be essentially larger than the renormalon picture predicts. Thus the asymptotic series for 
the pole mass seems to reach its critical order and starts to diverge. If this is a case then 
eq. ( ^41) gives a numerical estimate of a finite order sum of the asymptotic series. This value 
extracted from the sum rules can be used as an auxiliary parameter in the expressions for the 
physical quantities order by order of QCD perturbation theory. This, however, is sufficient 
for physical applications because the pole mass is not an observable. 

An independent NNLO analysis of T sum rules has been done in ref. where the Laplace 
transform of the polarization function was studied in the spirit of ref. while the explicit 
NNLO expression of the polarization function near the threshold has not been obtained. We 
found that for a given set of the parameters the numerical values of the theoretical moments 
obtained on the basis of our result for the polarization function is in a good agreement with 
the results of ref. . However, in ref. an essentially lower range of the soft normalization 
scale 1.5 GeV < /i < 3.5 GeV was used for the phenomenological analysis. This does not 
affect strongly the result for the strong coupling constant 0.96 < as{Mz) < 0.124 obtained 
in ref. but leads to larger value of the 6-quark pole mass 4.78 GeV < m^, < 4.98 GeV (the 
result of the constrained fit for the fixed of ref. 0) than the number given by eq. (p^). 
From Fig. 2 we, however, see that such a low normalization scale is unsuitable for a reliable 
estimate because the perturbation theory around the Coulomb solution diverges in this case. 

4 Inclusive B-meson semileptonic width and \Vcb\ ma- 
trix element. 

The theoretical expression for the inclusive semileptonic width of i?-meson up to the second 
order in the strong coupling constant and up to the first order of heavy quark expansion 
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reads 



r. = U {ml/ml) Cria.) 1 - - 2F, (ml/ml) ^ (25) 



where Fi{x) = 1 — 8x — 12x^lnx + — is the phase space factor, ^2(0;) = (1 — x^, fi^, 
and are the HQET parameters [p3| , |3^ . The hadronic matrix element of the gluon dipole 
operator /iq of the heavy quark is directly related to the masses of 5-mesons with different 
spin structure 

fil (1 + 0(l/m,)) = ^(M|. - M|) = 0.36 GeV\ 

The hadronic matrix element of the heavy quark kinetic operator /i^r suffers from larger 
uncertainty. For numerical estimates we use the result 

/i^ = (0.5 ±0.15) GeV^ 

obtained within QCD sum rules framework ||35|. The perturbative coefficient Cr(as) up to 



the second order in |T^, ^ is 



Cp = 1 - 1.67^^ - (8.4 ±OA)(^y (26) 



TT \ TT 



for the normalization point fi = ^Jrri\^c- 

Along with ordinary power corrections coming from heavy quark expansion an account 
for truncation of perturbative series for semileptonic width can generate further contributions 
to eq. (p5D which within the renormalon picture scale as AgcD/mfe. We discuss some details 
related to this issue later. 

One sees from formula (l25|) that the semileptonic width depends rather strongly on mf,. 
Therefore if is taken from a theoretical expression for some other experimental quantity 
it should be determined with a great accuracy for obtaining a reasonable precision of the 
\ych\ determination. In the previous Section we have described the moments of the vac- 
uum polarization function and have extracted the numerical value of the 6-quark pole mass 
entering the theoretical perturbation theory expressions for the moments. Being analyzed 
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independently, the perturbative series in as for the moments and for the width expressed 
in terms of the pole mass m;, seem not to enjoy a fast apparent convergence that can lead 
to a large uncertainty due to higher order contributions. In principle the apparent conver- 
gence of the finite order perturbation theory series can be changed by redefinition of the 
mass |38| or the coupling constant 0. However mb is not an observable and 

has no immediate physical meaning. Therefore it can be safely removed from relations be- 
tween physical observables. The idea of trading the unphysical parameters in favor of direct 
relations between observables constitutes now the most general trend in high precision phe- 
nomenological analyses, its application to the particular case of the 5-meson semileptonic 
decay width was discussed earlier (see e.g. p, ^ and references therein). In this Section we 
establish the direct relation between inclusive 5-meson semileptonic width and one of the 
moments obtained in the previous section. This relation is independent of any redefinition 
of the quark mass and of using different couplings so the accuracy of the relation reveals the 
actual precision of the approximation used for comparison of two physical observables. It 
happens that the convergence of the perturbation theory for this relation is rather fast and 
the subsequent approximations converge very well in a heuristic sence that the higher order 
approximation is close to the previous one. 

Our analysis consists in direct relating the factor m\ in eq. ( p5D to the rath moment of 
the T sum rules ^ 

where (dimensionful) moments M':^^ = Aifj^^/ml^ is a purely experimental quantity. The 
theoretical moment Ai^^ is a dimensionless quantity which depends on rrih only logarithmi- 
cally in a finite order in a^. Thus, the substitution of relation ( pTj) to eq. (^) substantially 
reduces its dependence on rrib though it also introduces an explicit uncertainty due to A^^^'^. 
We use eqs. to find for the mixing angle \Vcb\ 



( ]\4'exp\5/'in 

\V^b\ = (1927r^)V^ir, /-' > 



r (28) 
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where the functions 

K,, = (K'^)-^/- (^F, (ml/ml) Cr{a.) (l - - 2F, (ml/ml) ^) ' (29) 

accumulate theoretical information depending on rub, rric, as, etc. 

The function Fi(m^/m^) gives rather large theoretical uncertainty if masses of b- and c- 
quarks are considered as independent variables. However there is almost model independent 
constraint of the form 

mt-m,- III f-^ - + 0(l/mL) = - Md = 3.34 GeV (30) 

where Mb = 5.31 GeV, = 1.97 GeV denote the spin-average meson masses, e.g. Mb = 
\{Mb + 3Mb*). With this constraint the function Fi(m^/m^) becomes a function of a single 
variable Fi(mfc). Note that in such a setting the rrib dependence of the function Fi{mb) 
partly cancels the large dependence of the width. Furthermore if the relation (|27|) is 
used to express the 6-quark pole mass in the argument of the function Fi{mb) in terms of 
the moments of the spectral density then only logarithmic dependence on rrib appears in the 
right hand side of eq. (p^). 

Now we analyze the theoretical factor Kth numerically order by order in as- The result 
reads 

Kth\n=w = 1.366(1 + 0.088 + 0.028). (31) 

where the value of 6-quark pole mass in the fixed order in as with Coulomb resummation 
is found from eq. (pTj). For comparison, the perturbative series for rrib that follows from 
eq. (p7D and the series for Cr are 

mb\n=io = 4.71(1 - 0.001 + 0.021), (32) 

Cr = 1 - 0.146 - 0.064. (33) 

Thus we find that in the expansions of the theoretical moment (as well as rrib itself) and 
width expressed in terms of rrib the NNLO corrections are of the order of the NLO ones while 
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the perturbative series for the mixing angle, or the theoretical coefficient eq. (|3T|), converges 
much better. 

Let us discuss the problem of convergence in more details. Though there is no rigorous 
result on the asymptotic structure of the expansion (|26D it is widely believed (mainly due to 
consideration of the renormalon contribution) that this asymptotic series (as well as the series 
for the moments) starts to diverge already in some low orders of perturbation theory. The 
truncation of the asymptotic series at the optimal order results in the inherent uncertainty 
in these quantities which parametrically as large as KqcD/f^b within the renormalon picture. 
An argument in favor of this conjecture is the relatively large value of the NNLO corrections 
both to the moments and the semileptonic width. Since the detailed structure of the above 
asymptotic expansions remains unknown we cannot conclude that the divergence of the 
series for the moments and the series for the width cancels each other exactly. But because 
the terms in the perturbative expansion (^) decrease rapidly we can reasonably hope that 
the critical order where the series (|3ll) starts to diverge is higher than ones of eq. ( p2|) and 
eq. (^). Moreover, a partial cancellation of the divergences has been found within the 
renormalon picture [110, R^. 



As for numerics, we use the following central values for our experimental inputs (see |^ 
for more detail): 

BR(5 ^ XM) = 10.5%, tb = 1.55 ps, 

^ 3 ^ GeV-^ ^.(M^) = 0.118. 

With these numbers we obtain the value of the matrix element \Vcb\ 

, /BR(S ^ Vl-55ps\^ 
11^0.1 ^ 0.0423 0.105 j (^) 

x(l-0.01 °''-^^-)J^"« )(l±A„..)(l±A.) 

where /S.npt ~ 0.02 is the uncertainty in the nonperturbative contribution ||^ induced mainly 
by the uncertainty of /i^r in eq. (pO]). We have introduced also the uncertainty due to trun- 
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cation of the perturbative series for the Kth parameter Atr ~ 0.01 which is taken as a half 
of the last term in eq. (|3T|). 

The typical scale of uncertainty of key parameters is also indicated. Another important 
source of the uncertainty is the scale dependence of the theoretical moment and the exper- 
imental errors in the value of as{Mz) because of rather high sensitivity of the theoretical 
moment to In fact these uncertainties are closely related since as it has already been 
pointed out the scale dependence of A4^^{as{^), fi) is mainly due to the scale dependence of 
asifi). The pointed error bars roughly correspond to the interval 3.5 GeV < yU < 6.5 GeV 



at fixed as{Mz) = 0.118. The central value in eq. (34) does not change in the interval 
5 < n < 15 and we have chosen n = 10 for our final estimate by the reasons discussed in the 
previous section. 

The main part of the experimental uncertainty is related to the uncertainty in the exper- 
imentally measured inclusive semileptonic width. The experimental situation changes rather 
quickly and data are improving fast that means that the experimental uncertainties will be 
smaller (see e.g. 0]). The uncertainty in Mfg^ comes mainly from the uncertainties in 
leptonic widths of T resonances. It is about 5% and leads only to 0.3% uncertainty in 
\Vcb\ so we do not include it to the error bars in eq. (Q). 

Now the advantage of our approach becomes clear - large errors due to the uncertainty 
in rrib is now partly shifted to more direct experimental data. The expression for |Vc6| 
matrix element has a very weak dependence on so a possible uncertainty in 6-quark pole 
mass does not lead to any uncertainty in \Vcb\- Furthermore the terms in the perturbative 



expansion ( pTD decrease rapidly which is an indication that the higher orders corrections to 
the obtained result are small enough. 

Our result is in a good agreement with the previous estimate |V"c6| = 0.0419 f^. Our 
value, however, is somewhat larger than the estimate \Vcb\ = 0.039 of ref. @]. There is no 
much hope to reduce the uncertainty in the nonperturbative contribution. Thus, in our 
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opinion, the model independent result presented in the paper provides one with the most 
reliable and accurate estimate of the CKM matrix element | Vcb \ from the inclusive 5-meson 
semileptonic width. 

5 Conclusion. 

In this paper we have presented the determination of the 6-quark pole and MS mass and 
Cabibbo-Kobayashi-Maskawa matrix element |Vc6| from T sum rules and 5-meson semilep- 
tonic width in the NNLO in the strong coupling constant. In our opinion it seems to be 
a final result for these problems within the framework of analytical treatment: the next 
order approximation is too complicated to deal with analytically within QCD or NRQCD. 
In particular, the nonrelativistic approximation in next order in as should be supplemented 
by the real gluon radiation that takes the problem to the completely new level of complexity 
and makes it practically unsolvable within the framework presented here. 

The self-consistence of the T sum rules has been checked and the intervals for the relevant 
normalization scale and the moment numbers which can be used for a reliable estimates have 
been found. An asymptotic character of the perturbative expansion for the moments of the 
T system spectral density was discussed and a conjecture on the critical order where the 
series in for the moments starts to diverge has been made. We have also presented a 
new representation of the heavy quark vector current correlator near threshold based on the 
explicit formulas for the NNLO correction to the ^5*1 heavy quark bound state parameters^ 

We have constructed the direct relation between the moments of the T system spectral 
density and the inclusive S-meson semileptonic width up to the NNLO order. We have 
found that when the unphysical variable (the 6-quark pole mass) is removed the residual 
perturbation theory in the coupling constant for Tsi (or \Vcb\) works well and demonstrates 
^Whcn this paper was in preparation a paper appeared where the similar result was obtained. 
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nice convergence. 
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Appendix. 

A. The correction A^2^G due to the A2V part of the potential is (see eq. 



KAttJ Att 



\m=0 



+ (m + l)^i(m + 2) (Ci^ + 2L{k)Cl) + J(m)C, 



00 m— 1 

+2 ^ ^ F{m)F{n) ( -^1-^ (C^^ + 2L{k)Cl) + J(m, n)^^ 



m=ln=0 ^ "^-^ 



00 

+2 5] F(m) (Co' + L{k)Cl + (^^(A;) + K{m))Cl - {2^e + + 1)) (cl + 2L(fc)C2)) 



m=0 



+L(A;)Co' + ( -7EL(fc) + ]^L\k) ) + N{k)Cl 



where 



j e-^Li{z)Ll{z)\n\z)zdz 



I{m), m = n, 
J{m, n), m > n, 

J e-'L\{z) \n{z)dz = K{m) 



J X 

where dots stand for inessential k independent divergent part and 

/ vr^ 2 
/(m) = (m + 1) ^^(m + 2) - ^2(^ + 2) H 



^ 3 (m+l)2^ 

-2(^i(m + l)+7s), 

n+1 f , , 1 \ 

J(m, n) = 2 ^i(m - h 27^5 

m — n \ n + 1 J 



Tlx ~\~ 1 

+2 (*i(m - n + 1) - ^i(m + 1)), 

m — n 

K{m) = 2(^i(m + 1) + 7ij)^ + ^2{m + 1) - ^^(^ + 1) + 27I, 
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The correction A2^'*G due to the second iteration of AiV term 



\m=0 



{c'o + {'^{m + 2) + L{k))Cl 

oo m-1 „ _|_ 1 / / / 11 

/ oo m— 1 1 ^ I -| 



Vm=2 /=1 n=0 



(/ — n){m — n) 



oo m— 1 71—1 ; I 1 

B. The NNLO corrections to the square of the Coulomb ^Si heavy quark bound state wave 
function at the origin (eq. (|T5D) 



m 



3(m + l)2' 



^ -ex (f ^ . (M) (- . 



+7b + ^i("^ + 1) ^ 



(m + 1)^ 

A?V^ = (^)' (3(^0^ + L{m)Cl + L2(m)C2^) + (-1 - 27^ + ^ + vl/i(m + 2) 



-2(m + l)^2(m + \)){Cl + 2L(m)C2') + ( + 2ir(m) - 2^i(m + 2) 

V m + 1 



+2 E 7 ?7 rJ(m,n) + 2 ^ -J(n,m) C, 



n=0 



(n-m)(n + l) ' ^^^^'^i (n - m)(n + 1) 



2 I ' 



2 
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Z*^! (n + l)(m + l) /^i (.^. . ^, ,,1 n + l \ 
(n-m)^ V V 2n-mJ 

Cr m— 1 r) -I- 1 
n=0 

^ m + 1 ^, A ^,^1.2 ^ n + l ^ m + 1 



n=m+l 



(n-m)2 'I I 2 ^^{n-mY ^^^-^ {n - m) 



2 



l-^i n + l ^, , 1 ^ m+1 2 

2;^o(^-"^) 2„J^^Jn-m)3(n + 1) 

^^i^Y (/ + l)(n + l) (/ + l)(m + l) (/ + l)(m + l) 

n=l «=0 



(n — my{l — m)2 (n — my{l — m){n — I) (/ — my{n — m){n — I) 



^ - rn_l / + ^ (/ + 1)(^+1)2 

V {n-mf{l-mY {n - m)'^{l - m){n - l){n + 1) 
(/ + l)(m + l) \ g g (m + l)2 



+ 



(/ + l)(m + l)2 (m+l)2 



2 



where 



— my(l — m){n — l){n + 1) (/ — mY{n — m){n — I) 
n + l {n + lf 



U (m, n) = 3 + 



m + n + 2 — m) {n + m + 2) 



C The NNLO corrections to the Coulomb ^5*1 heavy quark bound state energy levels 
(eq. (0)) 

A F (Ca , 2 11 1 \ 



(m+1) \Cf 3 16(m+ l)y ' 



E„ = 2 (Co^ + L{m)Cl + L\m)Cl + ^i(m + 2)(C2 + 2L{m)Cl) 



^2 



(m+ 1) 

2 

47ry 



2 ' 



A^^E^ = + (L(m) - 2 + ^i(m + 2))Ci) (Cq^ + (L(m) + ^i(m + 2))C!; 
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+ 



-(7E + ^i(m + 2)) - 2^2(m + 1) - (m + 1)^3("^ + 1) (Ci 



1\2 



(m+l) 

The results for the heavy quark bound state parameters are obtained on the basis of correc- 
tions found in When this work was in preparation ref. |rD|| appeared where the correction 
to the ^5*3 heavy quark bound state wave function at the origin were computed up to NNLO. 
This calculation is consistent with our result for Aitp'^, Ad5ip'^, ANA,BF,A'^m A^2^iIj'^. 
The expression for A^^Vm 

was obtained in ref. |10[ in a different representation so direct 
comparison with our formula is difficult. At the same time numerically the results are in a 
good agreement. Our result for the NNLO correction to the heavy quark ^5*3 bound state 



energy levels coincides with the result of ref. [41 



Figure captions 

Fig. 1. The relative weight of the NLO corrections to the ground state energy 
(curve a) and to the square of the ground state wave function at the origin 
(|'?/'o(0)P)^^'^/(|'?/^o(0)P)'^'^ — 1 (curve b) as a function of the normalization point fi. 
Fig. 2. The relative weight of the NNLO corrections to the ground state energy 
j^NNLO I j^NLO _ ]^ (^Q^2;yg q,^ g^^d to the square of the ground state wave function at the origin 
(|'^o(0)n^^'^'^/(|'?/'o(0)n^'^'^ — 1 (curve h) as a function of the normalization point fi. 
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